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Abstract 

The  estimation  of  individual  instabilities  of  N  clocks,  compared  by  measuring  the  differences  of 
their  readings,  is  here  considered  without  assuming  a  priori  any  hypotheses  on  their  uncorrelation. 
Instabilities  of  the  N  docks  are  described  by  a  complete  ( non-diagonal )  NxN  covariance  matrix 
R.  Only  differences  of  dock  readings  are  available  in  order  to  estimate  R.  Statistical  processing  of 
these  data  allows  one  to  calculate  the  (N-l )x(N-l)  covariance  matrix  S  of  the  differences  relative 
to  the  N-th  (reference)  dock.  By  analyzing  the  relationships  tying  R  and  S,  several  pieces  of 
information  can  be  inferred  and,  in  particular,  the  conditions  for  the  validity  of  the  uncorrelation 
hypothesis  are  established.  The  estimation  of  R  from  S  is  not  unique:  in  any  case  R  must  be 
positive  definite.  A  theorem  states  that  R  is  positive  definite  if  and  only  if  its  determinant  is 
positive.  Nevertheless  infinitely  many  acceptable  choices  of  R  still  fulfil  the  condition  of  positive 
definiteness.  This  paper  shows  that,  by  increasing  the  number  N  of  compared  docks,  the  amount 
of  arbitrariness  in  estimating  R  is  reduced.  The  analysis  of  some  experimental  data  illustrates  the 
capability  of  the  method. 


1  INTRODUCTION 

The  evaluation  of  frequency  standard  instability  is  often  performed  by  comparing  N  clocks  and 
measuring  differences  in  their  signals.  Results  depend  on  the  simultaneous  contributions  of  ail 
N  standards  and  it  is  often  desirable  to  estimate  the  noise  contributions  of  any  single  units.  In 
the  past  years,  this  problem  has  been  considered  in  several  papers,  which  introduced  the  popular 
“3-cornered  hat”  method  [1],  successively  investigated  and  extended  to  N  clocks  [2,  3,  4,  5,  6]. 

The  application  of  these  methods  often  results  in  an  unpleasant  situation:  some  estimated  clock 
variances  turn  out  to  be  negative,  violating  the  positiveness  restriction  intrinsic  to  their  definition. 
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The  assumption  of  uncorrelation  of  clocks,  necessary  in  the  above  methods,  seems  to  be  too  strin¬ 
gent  and  not  always  supported  by  experimental  evidence  ([7],  see  also  [8]).  Moreover,  when  more 
than  3  clocks  are  compared,  the  uncorrelation  hypothesis  leads  to  the  formulation  of  an  overdeter¬ 
mined  linear  system  of  equations  [2,  3,  6,  9]  and  this  overdetermination  seems  incompatible  with 
the  inconsistency  of  some  results. 

Here  a  new,  consistent  and  more  general  mathematical  model  is  proposed,  suited  for  statistical 
processing  of  measured  differences.  The  entire  NxN  covariance  matrix  R  of  the  clocks  along  with 
the  (N-l)x(N-l)  covariance  matrix  S  of  the  differences  are  introduced.  Their  relationship,  and  the 
possibility  to  estimate  R  from  the  knowledge  of  S,  are  then  discussed.  A  diagonal  structure  of  R 
can  result,  a  posteriori,  from  the  most  desirable  situation  of  truly  uncorrelated  clocks.  This  model 
abides  by  the  CIPM  recommendation  [10]  in  which  it  is  recognized  that  “where  appropriate,  the 
estimated  covariances  should  be  given” . 

The  examination  of  matrix  S  helps  understanding  the  validity  of  the  uncorrelation  hypothesis. 
When  such  hypothesis  is  acceptable,  the  problem  can  be  reduced  to  the  classical  case  of  the  N- 
comered  hat,  while  in  the  other  cases,  the  complete  matrix  R  must  be  solved.  In  this  second  case 
the  problem  is  underdetermined  because  the  number  of  unknowns  exceeds  the  number  of  equations. 
In  any  case  there  is  an  important  constraint  which  binds  the  domain  of  acceptable  solutions:  the 
matrix  R,  as  any  covariance  matrix,  must  be  positive  definite.  A  theorem,  specifically  devised, 
is  here  presented  showing  that  R  is  positive  definite  if  and  only  if  its  determinant  is  positive. 
This  constraint  has  been  geometrically  interpreted  giving  an  insight  of  the  features  of  the  solution 
domain. 

The  interesting  case  of  an  increasing  number  of  compared  clocks  is  then  examined  showing  how  a 
larger  number  of  clocks  reduces  the  solution  domain  of  R. 

2  STATEMENT  OF  THE  PROBLEM 

Let  us  denote  by  xl  the  process  related  to  the  i-th  clock  and  xlm  its  realization  at  the  instant  tm. 
If  M  consecutive  equispaced  samples  of  the  process  are  considered,  they  can  be  arranged  in  the 
vector  x*  =  [xj,  x'MY,  where  superscript  T  denotes  the  transposition.  The  expected  value 

x1  of  process  x'  can  be  estimated  in  terms  of  the  elements  of  x*  as  xl  ~  (1  /M){x\  +  x2  +  ...  +  xlM) 
and  arranged  into  the  M-vector  x’  =  [5*,  5*, . . . ,  xl]T.  The  M  samples  of  the  N  clock  processes 
are  then  cast  in  the  M  x  N  matrix  X  =  [x1,  x2,  . . . ,  xw].  Similarly,  the  M  x  N  matrix  X  is 
introduced  as  X  =  [x1,  x2,  . . . ,  xw|.  With  these  definitions  the  covariance  matrix  of  the  processes 
xl,  i  =  1,  2, . . . ,  N  is  estimated  by  R: 

R=[l/(M-l)][X-Xf[X-X]  (1) 

The  diagonal  element  r«  (r*  >  0)  denotes  the  variance  of  the  i-th  clock,  while,  for  i  ^  j ,  r{}-  denotes 
the  covariance  between  the  i-th  and  the  j-th  clock.  In  order  to  obtain  a  good  estimation  of  the 
matrix  R,  it  is  required  that  the  number  of  samples  M  be  much  larger  than  the  number  N  of 
clocks. 

When  clocks  are  compared  by  measuring  differences  between  their  readings,  the  available  data  are 
the  time  differences  y'j  =  x1  —  xK  Also  yh  is  statistically  characterized  by  its  expected  value  y*J 
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and  by  its  variance:  the  estimates  of  these  values  are  tied  to  the  estimates  of  the  variances  r\;  and 
Vjj  of  the  two  single  clocks  and  of  their  covariance  ri3-. 

One  of  the  clocks,  for  instance  clock  #JV,  is  chosen  as  the  reference,  and  it  is  compared  at  M 
different  instants  with  clocks  #1,  #2,. . . ,  #(N-1)  giving  (JV  —  1)  distinct  time  difference  measure 
M-vectors  y1JV  =  x1  —  xw,  y2N  =  x2  —  xw, .  . . ,  N  =  x(JV_1)  —  x'v.  Under  the  assumptions  of 

contemporary  and  noiseless  measurements,  the  other  possible  measure  vectors  y'J  (i,  j  ^  JV )  are 
redundant,  since  they  can  be  obtained  as  linear  combination  of  those  involving  clock  #JV  and  they 
don’t  add  any  information,  as  proved  in  [11]. 

The  M  samples  of  the  IV  —  1  clock  differences  are  cast  in  the  M  x  ( JV  —  1 )  matrix  Y  =  [y1A\  y2N , . . . , 
yjv-i,  ,\j  Similarly  the  M  x  (jV  —  1)  matrix  Y  is  introduced  as  Y  =  [yiw,  yw, . . . ,  y^_1)’  w],  where 
yiN  is  the  M- vector  whose  elements  are  coincident  with  yiN. 

The  processes  y,v,  (i  =  1,  2, . . . ,  N  —  1)  are  statistically  characterized  by  the  estimation  of  the 
(N  —  1)  x  (JV  —  1)  covariance  matrix  S  defined  as: 


S  =  [1/(M-1)][Y-Y]t[Y-Y]  (2) 

A  generic  element  denotes  either  a  variance  (for  i  —  j)  and  it  is  always  positive  or  a  covariance 
(for  i  ^  j)  and  may  assume  any  real  value.  Note  that  the  index  N  of  the  reference  clock  has  been 
dropped  in  any  element  of  S. 

The  relationship  between  matrices  S  and  R  is  easily  determined  by  observing  that  Y  can  be 
derived  from  X,  according  to  the  following  relationship: 


Y  =  X  H 


(3) 


where  H  is  the  JV  x  (JV  —  1)  matrix: 


H  = 


u' 


(4) 


where  I  is  the  (JV  -  1)  x  (JV  —  1)  identity  matrix  and  u  is  the  (JV  -  l)-vector  [1,  1,. . . ,  1]T.  A 
similar  relationship  ties  the  corresponding  matrices  X  and  Y. 

From  (1),  (2),  and  (3),  the  covariance  matrix  S  can  be  expressed  in  terms  of  R  as: 


s  =  htrh 


(5) 


For  sake  of  generality  here  we  discuss  the  covariance  matrices  of  the  processes  themselves,  but  the 
same  properties  hold  also  for  any  covariance  matrix  defined  from  prefiltered  data  as  is  generally 
the  case  in  the  estimation  of  clock  instability,  where  the  Allan  variance  is  used. 

3  CONSIDERATIONS  ON  R  AND  S 

By  taking  into  account  the  symmetry  of  R  and  S  it  appears  that  R  is  defined  by  JV(JV  -  1)  scalars 
and  S  by  (JV  —  l)JV/2  ones,  then  the  knowledge  of  S  is  not  sufficient  to  fix  a  unique  estimation 
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of  R  because  oo'v  solutions  satisfy  (5)  (underdetermined  problem).  A  way  to  isolate  the  N  free 
parameters  consists  in  partitioning  R  according  to: 


(6) 


where  R  is  the  leading  (N—  1)  x  (N- 1)  submatrix;  r  =  [riw,  ,  rN-h  w]3"  is  the  (N  —  l)-vector, 

grouping  the  covariances  involving  the  iV-th  clock  and  rNN  is  the  variance  of  the  JV-th  clock.  With 
this  partition,  and  by  substituting  (4),  equation  (5)  can  be  rewritten  as: 


s  =  [!  -  u] 


=  R  +  rjW[u  u]r  —  u  rr  —  r  uT 


This  equation  shows  that  R  can  be  uniquely  reconstructed  from  S  if  the  (co)variances  related 
to  the  iV-th  clock,  i.e.  r  and  rNN,  are  known.  In  fact,  (7)  can  be  transformed  in  the  following 


expression: 

R  =  S  -  rWJV[u  ur]  +  u  rr  +  r  uT 


(8) 


With  a  deeper  look  into  (7),  some  qualitative  information  about  R  can  be  immediately  deduced 
when  S  is  known: 


1.  If  the  N  clocks  are  uncorrelated,  S  assumes  the  form: 

Ul+T/VAT  rNN  rNN 

Tnn  ^22  +T 'nn  •••  rNN 

Tun  Tjvj v  Ty-i,  w-i  +  T"nn 

All  the  covariances  (with  i  ^  j)  are  equal  and  positive,  because  ~  rNN  >  0.  Moreover 
each  variance  >  s{j  (with  i  ^  j).  These  conditions  are  necessary  to  validate  the  assumption 
of  uncorrelation  of  the  N  clocks.  Unfortunately  they  are  not  sufficient  to  conclude  that  the 
clock  are  uncorrelated  because,  together  with  a  unique  diagonal  solution,  infinitely  many 
other  non-diagonal  matrices  R  would  drive  to  the  same  matrix  S,  but  if  these  conditions  are 
satisfied  one  can  “reasonably”  assume  the  uncorr elation. 

2.  If  the  reference  clock  is  “quasi-ideal”  (that  is  rNN  <§;  rti  for  i  =  1,  2, . . . ,  N  —  1)  and  all  the 
clocks  are  uncorrelated,  then  -C  su  for  any  i  ^  j  and  S,  as  well  as  R,  can  be  considered 
diagonal.  So,  if  S  is  almost  diagonal,  the  reference  clock  is  of  high  quality.  In  this  case,  the 
submatrix  R  almost  coincides  with  S  and  r  =  0,  rNN  =  0. 

3.  For  the  same  reason,  when  the  hypothesis  of  uncorrelation  holds,  if  the  reference  clock  is 
changed  and  S  is  computed  again,  it  gives  an  idea  of  which  one  of  the  clocks  is  less  noisy, 
because  it  results  in  a  matrix  S  with  minimum  off-diagonal  terms. 

4.  If  any  Sij  (with  i  ^  j )  is  negative,  then  the  uncorrelation  hypothesis  is  certainly  to  be  excluded 
at  least  between  the  (i,  IV)-th  or  the  ( j ,  iV)-th  pair  of  clocks. 
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5.  If  all  the  terms  (with  i  ^  j)  are  equal  and  positive  except  one,  for  instance  sH(with 
k  ^  l),  which  is  positive  but  differs  from  the  others,  that  indicates  that  all  the  clocks  can  be 
considered  uncorrelated  except  the  pair  (k,  l). 

6.  If  the  terms  (with  i  ^  j  ^  k)  are  coincident  and  positive,  and  the  terms  (with 
i  ■£.  j,  i  —  k,  or  j  =  k)  are  coincident  but  different  from  the  previous  ones,  the  fc-th  clock  is 
possibly  correlated  with  the  reference  one  while  all  the  others  are  uncorrelated. 

From  the  theoretical  point  of  view,  the  above  considerations  suggest  clear  interpretations.  By 
handling  experimental  data,  it  can  be  difficult  to  recognize  the  points  above  because  of  round-off 
error,  not  perfect  contemporaneity  of  measurement,  low  noise  introduced  by  the  measurement 
system  but,  above  all,  the  low  number  of  statistical  samples  that  give  a  considerable  uncertainty 
bar  on  the  estimates. 

4  POSITIVE  DEFINITENESS  OF  S  AND  R 

The  estimation  of  the  complete  matrix  R  asks  for  a  suitable  choice  on  the  A  free  parameters  in 
(8).  In  any  case,  there  is  an  important  constraint  which  bounds  the  solution  domain  V  in  the 
space  of  the  A  free  parameters  and  which  guarantees  a  significant  result:  as  any  covariance  matrix 
[13],  the  estimated  covariance  matrix  R  must  be  positive  definite. 

Some  significant  properties  of  positive  definite  covariance  matrices  are  here  reminded.  Let  A  be 
an  A  x  A  symmetrical  matrix  and  Aw  (n  =  1,  2, . . . ,  A)  be  the  leading  n  x  n  submatrices  of  A, 
extracted  from  A  by  considering  the  elements  belonging  to  the  first  n  rows  and  n  columns  (with 
Aw  =  A).  The  matrix  A  is  positive  definite  if  and  only  if  [13]: 

| A”) |  >  0  for  n  =  1,  2,...,  A  (10) 

where  |  .  |  denotes  the  determinant. 

Then  the  positive  definiteness  of  A  requires  the  validation  of  A  scalar  inequalities.  In  the  case  of 
R  and  S,  since  they  are  both  positive  definite  and  they  are  linked  by  (5),  a  theorem  holds  [llj: 

Theorem  1:  The  A  x  Asymmetrical  matrix  R,  related  to  the  known  positive  definite  ( A  —  1)  x 
(A  —  1)  symmetrical  matrix  S  by  the  relationship  (5),  is  also  positive  definite  if  and  only  if  [R|  >  0. 

This  theorem  allows  one  to  verify  the  positive  definiteness  of  R  by  considering  only  one  of  the 
inequalities  (10);  moreover  this  inequality  can  be  geometrically  interpreted.  In  fact,  the  relationship 
|R|  =  0  appears  to  be  an  equation  of  second  degree  in  the  A  —  1  elements  of  r  and  in  rNN. 

With  the  partition  of  R  used  in  (6),  it  can  be  demonstrated  [11]  that  its  determinant  can  be 
expressed  as: 

|Rj  =  ISKrjv*  -  [r  -  r^uf  S_1[r  -  r**u])  (11) 

with  u  defined  in  (4) .  The  boundary  of  the  solution  domain  T>  is  therefore  described  by  the  surface 
|R|  =  0: 

[r  -  rwwu]rS  x[r  -  rwu]  =  rNN  (12) 

If  rNN  is  fixed,  this  expression  represents  an  ellipsoid  in  the  (A  —  1) dimensional  space  described  by 
the  variables  r1A,,  r2N, . . . ,  rw_u  N  because  the  matrix  S_l  of  the  quadratic  form  is  positive  definite. 
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This  ellipsoid  is  centered  at  the  point  of  coordinates  (r,w  =  rNN,  r2N  =  rNN,... ,  rN-x<  N  =  rNN), 
the  principal  axes  are  rotated  and  the  squared  length  of  the  axes  is  proportional  to  rNN.  As  far 
as  the  parameter  rNN  is  concerned,  the  domain  T>  takes  the  form  of  an  elliptic  hyper-paraboloid 
in  the  A- dimensional  space.  The  intersection  with  a  plane  (rNN  —  constant)  leads  to  an  (A  —  1)- 
dimensional  ellipsoid,  the  increase  of  rNN  increases  the  surface  and  shifts  the  center  of  the  ellipsoid 
(12),  and  it  doesn’t  rotate  its  axes.  Fig.  1  shows  an  example  of  the  domain  T>  in  the  case  of  A  —  3. 

The  constraint  of  positive  definiteness  of  R  bounds  the  domain  of  choice  of  the  A  free  parameters  r 
and  rNN  and  only  the  set  of  values  of  riw,  r2,v, . . . ,  rN_lt  N,  rNN  situated  inside  the  surface  described 
by  (12)  guarantees  acceptable  solutions. 

5  EFFECT  OF  AN  INCREASING  NUMBER  OF  COMPARED 
CLOCKS 

Knowing  the  expression  of  the  boundary  of  the  domain  T>,  it  can  be  interesting  to  test  how  it  is 
modified  by  varying  the  number  of  clocks  involved  in  the  measurements.  Within  a  set  of  A  clocks 
let’s  individuate  the  A'-th  clock  as  the  reference  one  and  two  subset  Ca  of  A“  clocks  and  Cb  of  Nb 
clocks.  C  denotes  the  union  of  Ca  and  Cb  and  A  =  A“  +  Nb  +  1.  With  the  above  notation,  the 
A  x  A  covariance  matrix  R  of  the  A  clocks  is  partitioned  as  follows: 


R  = 


Rm  R°*  ra 
[R"*f  R“  r* 

[r*f  [r‘r  r„ 


(13) 


where  Raa  and  R66  are  the  covariance  matrices  related  to  the  clocks  in  Ca  and  Cb  respectively. 
Ra'>  contains  the  covariances  between  clocks  in  Ca  and  Cb  versus  the  A'-th  clock;  they  represent  a 
partition  of  the  vector  r  introduced  in  (6).  The  element  rNN  is  the  variance  of  the  A-th  clock. 


Analogously,  the  (A  —  1)  x  (A  —  1)  covariance  matrix  S  of  the  differences  can  be  partitioned  as: 


S  = 


S'“ 

[S*]T 


(14) 


where  S“  and  S'*  are  the  covariance  matrices  of  the  differences  related  to  the  clocks  in  Ca  and  Cb 
respectively. 

If  the  clocks  in  Ca  and  C  (jointly  with  the  reference  one)  are  considered,  the  respective  solution 
domains  T>a  and  T>  are: 

[ra  -  rN,vu°]T[S“]  ‘1[ra  -  rWNua]  <  rNN  (15) 

[r  -  rJvwu]:rS-l[r  -  rwwu]  <  rNN  (16) 

where  u“  is  the  Na- vector  [1,  1, . . . ,  \]T . 


Expression  (16)  can  be  rewritten  in  terms  of  the  submatrices  of  S  in  (14).  In  fact,  the  vector 
r  —  rjVlVu  can  be  partitioned  as: 


[r  -  ^u] 


r°  -  rNNua 
rb  -  t-jvjvU6 


(17) 
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where  is  the  iV^-vector  [1,  1, . . . ,  l]r  and  the  inverse  of  S  can  be  written  as: 


s-1  = 


[goaj-i  +  [gocj-1  g ai  [S“j_1 

-  [D“]_1  [S'36]7"  [S“]_1 


-  [S0*]-1  S”4  [D44]”1 
[D44]-1 


where  the  matrix  D44  is  defined  as  [12]: 

D44  =  S44  -  [S’*]T[S“]-1Sa* 


(18) 


(19) 


Now  the  idea  is  to  compare  the  extension  of  the  domains  delimited  by  the  elliptic  paraboloids  IP" 
and  T> ,  defined  by  (15)  and  (16),  respectively.  The  domain  is  defined  in  the  space  of  Na  -f  1  vari¬ 
ables  r1N,  r2N, . . . ,  rN«N,  rNNi  while  V  is  defined  in  the  space  of  N  variables  r1N,  r2N, ... ,  rN_1>lV,  rM. 
To  the  aim  of  comparison,  it  is  necessary  to  individuate  a  sub-domain  V  obtained  from  T>  by  a 
suitable  projection  in  the  subspace  of  Na  +  1  variables  r1JV,  r2jV,  •  •  • ,  tvn,  rNN. 

Domain  T>  is  obtained  by  the  projection  of  the  intersection  of  T>  with  the  A^-dimensional  subspace 
defined  by  the  hyperplane  tvjVu&  =  rb.  Since  this  choice  corresponds  to  the  largest  intersection  of 
V,  V  results  in  the  largest  possible  subdomain  in  the  (Na  +  l)-dimensional  subspace  and  therefore 
it  yields  the  least  favourable  case  in  the  comparison  of  the  two  domains. 

By  substituting  (17),  (18),  (19)  and  =  rf’  in  (16),  the  equation  defining  domain  f>  is  obtained 

as: 


[ra  - r^uTtSTV -  rVlVu“]  +  [ra  - rlWuan[S“] 


Ml-lfcoblT 


[S“]  ^[r13  —  r^ua]  —  Tun  (20) 


By  comparing  V  (20)  with  Va  (15),  it  appears  that  the  only  difference  is  in  the  second  quadratic 
form  in  the  l.h.s.  of  (20).  From  the  properties  of  the  definite  positive  matrices,  since  S  and  its 
submatrix  S““  are  positive  definite,  so  are  their  inverses,  and  it  can  be  demonstrated  that  also  the 
matrix  [S“a]_1S“4[D4*]-I[S',4]r[Sa“]_1  of  the  second  quadratic  form  in  (20)  is  positive  definite.  So  it 
necessarily  results  that  T>  is  contained  inside  T>a.  As  a  result,  an  increased  number  of  compared 
clocks  reduces  the  solution  domain  and  therefore  the  amount  of  arbitrariness  in  the  determination 
of  the  covariance  matrix  R,  even  in  the  worst  considered  case,  i.  e.  with  the  largest  possible 
intersection.  It  must  be  stressed  that  the  reduction  ratio  of  the  solution  domains  depends  only  on 
the  elements  of  matrix  S  and  not  on  the  value  of  rNN .  Intuitively,  if  the  clocks  in  Cb  are  completely 
correlated  with  any  clock  in  Ca  or  with  the  reference  one,  their  measurement  doesn’t  add  any 
further  information  to  the  problem  and  the  reduction  of  the  solution  domain  will  be  very  poor.  On 
the  contrary,  the  addition  of  independent  clocks  can  result  in  a  significant  shrinking  of  the  solution 
domain. 


6  EXAMPLE 

Five  commercial  cesium  clocks  have  been  compared.  Three  of  the  clocks  are  maintained  at  IEN, 
Turin,  Italy  (Cs4,  Cs5,  Cs6),  the  other  two  clocks  are  maintained  at  ISPT,  Rome,  Italy  and 
Telespazio,  Matera,  Italy  and  data  are  regularly  transmitted  to  the  IEN  for  traceability.  The 
comparison  at  distance  is  obtained  with  GPS  or  TV  link.  An  appropriate  smoothing  procedure 
has  been  applied  to  reduce  the  additive  noise  due  to  the  synchronization  link  and,  for  integration 
times  larger  than  five  days,  this  noise  can  be  assumed  to  be  negligible.  The  Allan  covariance  matrix 
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has  been  computed  for  different  integration  times,  with  overlapping  samples,  by  using  one  year 
(1991)  of  daily  measurement  data.  (The  Allan  covariance  matrix  is  an  extension  of  the  well-known 
Allan  variance  in  the  sense  of  (1)).  In  the  present  example  the  set  Ca  is  composed  of  clocks  Cs4 
and  Cs5  of  IEN,  the  set  Cb  is  composed  of  the  two  external  clocks  and  the  clock  Cs6  of  IEN  is 
the  reference  one.  When  the  clocks  in  Ca  jointly  to  the  reference  are  used,  the  2x2  matrix  S““  is 
obtained  while,  with  the  five  clocks,  the  4x4  matrix  S  is  obtained. 

To  investigate  how  the  increased  number  of  clocks  reduces  the  solution  domain,  V  is  compared 
with  T>°  in  the  plane  (r15,  r2S) ,  with  a  fixed  value  for  the  variance  r65.  Here  we  report  only  some  of 
the  most  significant  examined  cases. 

The  first  case  regards  the  Allan  covariance  matrix  evaluated  for  integration  time  r  =  10  days.  The 
following  matrices  S“  and  S  are  obtained,  where  numerical  values  have  been  scaled  by  105S: 


11.8 

4.78 

8.64 

7.10 

4.78  ' 

c _ 

4.78 

9.34 

7.38 

1.31 

9.34 

o  — 

8.64 

7.38 

96.7 

-23.7 

7.10 

1.31 

-23.7 

218 

(21) 


Since  matrix  S“  satisfies  conditions  (9),  the  uncorrelation  hypothesis  is  acceptable  for  the  clocks 
in  Ca  and  the  reference  one.  We  can  choose  rK  =  s12,  ris  =  0,  and  rM  =  0  driving  to  the  classical 
solution.  By  choosing  =  s12,  the  two  domains  Va  and  T>  delimiting  the  possible  choices  for 
r15  and  t’25  are  compared  in  Fig.  2.  As  expected,  the  classical  solution  corresponding  to  rl5  =  0 
and  r25  =  0  is  inside  the  domain  P,  and  it  is  also  inside  the  domain  T>.  It  also  means  that  from 
the  comparisons  of  the  five  clocks  it  results  that  the  clocks  in  Ca  and  the  reference  one  may  be 
considered  uncorrelated.  That  is  not  the  case  for  the  other  clocks  because  the  complete  matrix  S 
has  not  all  equal  and  positive  off-diagonal  elements. 

From  Fig.2  it  can  be  verified  that  the  addition  of  the  (co)variances  of  the  clocks  in  Cb  slightly  shrinks 
the  solution  domain  and  therefore  also  the  arbitrariness  in  the  choices  of  the  free  parameters 
and  fa  is  reduced. 

The  second  case  concerns  the  instability  estimation  of  the  same  clocks  for  r  =  100  days.  The 
following  matrices  S“  and  S  are  obtained: 


S00  = 


104  2d.3 
20.3  16.1 


104 

20.3 

97.4 

103 

20.3 

16.1 

19.8 

-41.6 

97.4 

19.8 

97.3 

88.9 

103 

-41.6 

88.9 

433 

(22) 


where  the  values  are  again  scaled  by  1028. 

Here  the  conditions  (9)  are  not  satisfied  neither  in  the  comparison  of  three  clocks  (SM),  nor  with 
five  clocks  (S).  The  solution  domain  of  the  covariances  ris  and  rM  is  represented  in  Fig.  3.  Here 
we  see  a  significant  reduction  of  V  with  respect  to  T>a  when  five  clocks  are  compared.  The  clocks 
in  Ca  add  so  much  information  that  the  range  of  acceptable  values  for  r15  and  r25  is  reduced  by  a 
factor  of  5.  Inside  the  domain  P,  the  arbitrariness  in  choosing  ru  and  is  so  much  reduced  that 
it  seems  not  so  much  crucial  the  criterion  to  be  used  for  the  final  choice.  This  situation  is  also 
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depicted  in  Fig.  4  where  the  two  3-dimensional  elliptic  paraboloids  describe  the  boundaries  of  the 
solution  domains  for  r15,  r^,  and  r5i. 

7  CONCLUSIONS 


This  work  has  considered  the  estimation  of  clock  instabilities  when  3  or  more  clocks  are  compared 
by  measuring  Signal  differences.  From  the  measure  covariance  matrix  S  several  information  about 
the  instability  of  the  individual  clocks  can  be  inferred  and,  in  particular,  the  uncorrelation  hypoth¬ 
esis  can  be  accepted  or  rejected.  The  estimation  of  the  individual  clock  covariance  matrix  R  from 
the  knowledge  of  S  is  an  underdetermined  problem  and  an  appropriate  optimization  criterion  ought 
to  be  formulated  to  fix  the  remaining  free  parameters,  but  the  solution  domain  is  constrained  by 
the  request  of  positive  definiteness  for  the  resulting  R.  Such  a  constraint  has  been  deeply  analyzed 
and,  with  a  theorem,  it  has  been  possible  to  give  a  simply  geometrical  illustration.  Moreover  the 
importance  of  the  comparison  of  a  larger  number  of  clocks  is  outlined,  showing  how,  in  some  cases, 
it  results  in  a  significant  reduction  of  the  arbitrariness  in  estimating  clock  instability. 
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Fig.l:  Solution  domain  D  corresponding  to  a  particular  matrix  S  obtained  from  the  comparison  of 
3  different  clocks. 

(a)  the  domain  D  in  the  space  of  variables  r13,  r23,  and  r33. 

(b)  boundaries  of  D  for  some  fixed  values  of  r33. 


r25 


Fig-2:  Domains  If  (solid  line)  and  D  (dashed  line)  in  the  plane  (r15,r15)  for  integration  time  x  =  10 
days  and  r55  -  4.78*  10'28. 


QUESTIONS  AND  ANSWERS 


Question;  Is  it  possible  for  the  domain  to  shrink  to  zero? 

P.  Tkvella,  Institute  Elettrotecnico  Nazionale:  Yes,  in  fact  for  example,  I  suppose  that  one 
correlation  is  possible.  You  just  shrink  to  the  region  when  all  the  covariances  term  are  zero. 
You  just  have  no  more  degrees  of  freedom. 
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